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ABSTRACT 

We investigate the impact of high optical depth on the HI saturation observed in the Perseus molec¬ 
ular cloud by using Arecibo HI emission and absorption measurements toward 26 radio continuum 
sources. The spin temperature and optical depth of individual HI components are derived along each 
line-of-sight, enabling us to estimate the correction for high optical depth. We examine two differ¬ 
ent methods for the correction, Gaussian decomposition and isothermal methods, and find that they 
are consistent (maximum correction factor ^ 1.2) likely due to the relatively low optical depth and 
insignificant contribution from the diffuse radio continuum emission for Perseus. We apply the cor¬ 
rection to the optically thin HI column density on a pixel-by-pixel basis, and find that the total HI 
mass increases by ~10%. Using the corrected HI column density image and far-infrared data from the 
IRIS Survey, we then derive the H 2 column density on ^0.4 pc scales. For five dark and star-forming 
sub-regions, the HI surface density is uniform with Ehi 7-9 Mq pc“^, in agreement with the min¬ 
imum HI surface density required for shielding H 2 against photodissociation. As a result, Eh 2 /Shi 
and Ehi + Eh 2 show a tight relation. Our results are consistent with predictions for H 2 formation in 
steady state and chemical equilibrium, and suggest that H 2 formation is mainly responsible for the 
Ehi saturation in Perseus. We also compare the optically thick HI with the observed “GO-dark” gas, 
and find that the optically thick HI only accounts for ^20% of the “CO-dark” gas in Perseus. 

Subject headings: ISM: individual objects (Perseus) — ISM: molecules — radio lines: ISM 


1. INTRODUCTION 

As the most abundant molecular species in the uni¬ 
verse, molecular hydrogen (H 2 ) is the main constituent of 
giant molecular clouds, th e exclusive birthplaces of stars 
(e.g., iKennicutt fc Evaii^ 2012). Recent observations 
of galaxies at both low and high redshifts have shown 
that star formation rates are strongly correlated with 
H 2 surface densities (Eh 2 ), tbe relation g enerally k nown 
as the “K ennic utt-Schmidt law” l e.g.. iSchmidtl 1959; 
Kennicutfl 198 9: iBigiel et al.l 2008: IWilson et al. 2009; 


Tacconi et al.l 2010: ISchruba et al] 2011: iGenzel et ^ 


2013). This suggests that physical processes responsible 
for the atomic-to-molecular hydrogen (HI-to-H 2 ) transi¬ 
tion play a key role in the evolution of galaxies. 

Observationally, the HI-to-H 2 transition has been stud¬ 
ied via ultraviolet (UV) absorption measurements along 
many random line s of sight throu gh the Galaxy (e.g., 

I Savage et al.l 1977: iRachford et al.l 2002: iCTlmon et ^ 
2006). For these measurements, either early-type stars or 
active galactic nuclei were used as background sources, 
and HI and H 2 column densities, A^(HI) and Af(H 2 ), 
were estimated from Lyman-alpha and Lyman-Werner 
(LW) band absorption. The UV studies probed the 
II 2 mass fraction /h 2 = 2 V(H 2 )/[V(HI) -|- 2 V(H 2 )] 
ranging from ^10“® to ^10“^, and found that /h 2 
sharply increases at the total gas column density iV(H) 
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= V(HI) -b 2V(H2) of ~(3-5) X lO^o cm-^. Ad¬ 
ditionally, the HI-to-H 2 transition has been indirectly 
inferred from the flattening of the relation between 
the HI column density and a tracer of total gas col¬ 
umn dens ity (e.g., far- infrare d (FIR) or hydroxide (OH) 
emissi on: iReach et al.l 19 94: iMeverdierks fc HeithausCT 
1996: iDouglas fc TavloJ 2007: iBarriault et al.l 2010; 
iLisztl 2014). These studies found that the HI col¬ 
umn density saturates to ^(5-10) x 10^° cm“^, sug¬ 
gesting the presence of H 2 . The HI saturation has 
also been fo und in extragal actic o bservations on ^kpc 
scales (e.g., I Wong fc Blitz! 2002: iBlitz fc Rosolowskvl 

2006: lLerov et al.l 2008: lWong et all 20091. 

Theoretically, the HI-to-H 2 transition has been in¬ 
vestigated as a central process in photodissociation 
regions (PDRs). In PDRs, the interstellar medium 
(ISM) is predominatly atomic, and the molecular gas 
is only found in well-shielded regions where dissociat¬ 
ing UV photons are sufficiently attenuated. Many stud¬ 
ies have been presented with different treat ments of 
chem istry, geometry, and radi ative transfer le.g.. iSoitzed 
1948:lGoulcl fc Saloeteil 196 3: iGlassgold fc Langed 1974; 
Ivan Dishoeck fc BlackI 19 86: ISternberg 1988; ElmegreenI 
1993; Draine fc Bertoldil 1996: ^aans fc Neufeldl 1997: 


Browning et al. 2003: lGoldsmith et al.l 2007~ Liszt! 2007: 
Krumholz et alJ 2009 : iGlover et al.l 2010: lOffner et al.l 


2013: ISternberg et al.l 2014), and an ex cellent summary 
of the se studies was recently provided bv ISternberg et all 

(IMl . 

Among the many studies, the iKrumholz et al.l (|2009r) 
(KMT09 hereafter) model has recently been tested with 
a variety of G alactic and extrag alacti c observation s (e.g., 
Bolatto et al.l 201 1: iLee et al. 2012; IWeltv et al.l 2012; 


Wong et al.l 2013: lMotte et al.l 2014), thanks to its sim- 
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pie analytic predictions that allow a comparison with di¬ 
rect observables and an extrapolation of the model over a 
wide range of ISM environments. In the KMT09 model, a 
spherical cloud is illuminated by a uniform and isotropic 
radiation field, and the H 2 abundance is computed based 
on the balance between the rate of formation on dust 
grains and the rate of dissociation by UV photons (chem¬ 
ical equilibrium). The authors derived two dimensionless 
parameters that determine the location of the HI-to-H 2 
transition in the cloud, resulting in the following impor¬ 
tant predictions. First, they found that H 2 formation 
requires a certain amount of HI surface density (Shi) for 
shielding against the dissociating radiation field. Inter¬ 
estingly, this shielding surface density primarily depends 
on metallicity, and is expected to be ~10 Mq pc“^ (cor¬ 
responding to A^(HI) ~ 1.3 X 10^^ cm“^) for solar metal¬ 
licity. Second, the H 2 -to-HI ratio, i?H 2 = Sh 2 /Shi, was 
predicted to linearly increase with the total gas surface 
density. This is because once the minimum HI surface 
density is obtained for shielding H 2 against photodissoci¬ 
ation, all additional hydrogen is fully converted into H 2 
and the HI surface density remains constant. As a result, 
i?H2 is simply a function of metallicity and total gas sur¬ 
face density. Another interesting feature of the KMT09 
model is that the ISM is “self-regulated” in that pres¬ 
sure balance between the cold neutral medium (CNM) 
and the warm neutral medium (WNM) determines the 
ratio of the UV intensity to the HI density. 

Aiming at testing the KMT09 model on sub-pc scales, 
we have re cently focused on the Perseus molecular cloud 
(|Lee et al.l 2012). Perseus is one of the nearby molecular 
clouds in th e Gould’s Belt, a nd is located at a distance 
of ~300 pc (|Herbig fc Jones! 1983; ICernisI 1990). It has 
a projected angular size of ~6° x 3° on the sky (based 
on the CO emission)®, and lies at high Galactic latitude 
b ~ —20°, resulting in relatively simple HI spectra com¬ 
pared to other molecular clouds in the Galactic plane. 
With a total mass of ^2 x 10^ Mq (ISancisi et al.l 1974; 
iLada et ^ 2010), Perseus is considered as a low-mass 
molecul ar cloud wit h an intermediate level of star for¬ 
mation (iBallv et al.l 2008). To test the KMT09 model, 
we derived Shi and Eh 2 images using HI data from the 
Gal actic Arecibo L-ban d Feed Array HI S urvey (GALFA- 
HI: IStanimirovic et al.l 2006: IPeek et al.l 2011) and FIR 
data from the Improved Reprocessing o f the IRAS Sur¬ 
vey HRIS: lMiville-Deschenes fc Laeachi 2005b The final 
images were at ^0.4 pc resolution, and covered the far 
outskirts of the cloud as well as the main body. We found 
that the HI surface density is relatively uniform with Ehi 
^ 6-8 M 0 pc“^ for five dark and star-forming regions in 
Perseus (B5, BIE, Bl, IG348, and NGC1333). In ad¬ 
dition, the relation between i?H 2 and Ehi + Eh 2 on a 
log-linear scale was remarkably consistent for all individ¬ 
ual regions, having a steep rise of Rh 2 at small Ehi-)-Eh 2 , 
a turnover at i?H 2 ~ 1, and a slow increase toward larger 
i?H2- All these results were in excellent agreement with 
the KMT09 predictions for solar metallicity®, suggesting 
that the KMT09 model captures well the fundamental 
physics of H 2 formation on sub-pc scales. 

® In this paper, ^^CO(J = 1 —I 0) is quoted as CO. 

® Perseus has solar m etallicity llCo nzalez Hernandez et all 
2009). See Section 7.2.1 of ILee et al.l II2012I1 for a detailed dis- 
cussion. 


The observed HI saturation in Perseus, however, could 
alternatively result from the high optical depth HI. When 
the HI emission is optically thick, the brightness temper¬ 
ature (Tb) becomes comparable to the kinetic tempera¬ 
ture (Tk). As a result, the HI surface density saturates 
in the op tically thin approximation (which we used in 
iLee et all 2012) since Ehi oc Tb ~ Tk, and is underes¬ 
timated. As the constant HI surface density is the key 
prediction from KMT09, it is critical to evaluate how 
much of the HI column density distribution is affected 
by the optically thick HL 

In this paper, we assess the impact of high opti¬ 
cal depth on the observed HI saturation in Perseus by 
using HI emission and absorption measurements ob- 
tained toward 26 b ackground radio continuum sources 
(|Stanimirovic et HI 2014; Paper I hereafter). These ob¬ 
servations provide the most direct way to measure the 
high optical depth HI, allowing us to derive the “true” 
total HI column density distribution. Specifically, we 
use the HI emission and absorption spectra to estimate 
the correcton factor for high optical depth, and apply 
the correction to the HI column density image computed 
by ILee et all (j2012t) in the optically thin approxima¬ 
tion. We take this localized approach rather than using 
data from existing all-sky surveys le.g.. lHeiles fc Trolandl 
2003a; HT03a hereafter) in order to treat all spectra 
uniformly for the velocity range of Perseus and con¬ 
sider the possibility that GNM/WNM properties may 
vary with ISM environments (e.g., metallicity, star for¬ 
mation r ate, etc.) as expec ted fr om theoretical mod- 
els le. g.. iMcKee fc Ostrikeil 1977: iKovama fc In^utsuke] 
2002: IWolfire et al.l 200“^ lAudit fc Hennebelld 2005: 
iMac Low et al.l 2005: iKim et al.l 2013L 

This paper is organized in the following way. We start 
with a summary of previous studies where various meth¬ 
ods have been employed to derive the correction for high 
optical depth (Section [5|). We then provide a description 
of the data used in this study (Section |31) . In Section 
we estimate the correction factor for high optical depth 
using two different methods, and compare our results 
with previous studies. In Sections [S] and [H we apply the 
correc tion to the HI column density image from iLee et al.l 
(I2OI2II on a pixel-by-pixel basis, and revisit the HI satu¬ 
ration issue by rederiving the H 2 column density image 
and comparing our results with the KMT09 predictions. 
We then investigate whether or not the optically thick HI 
can explain the observed “GO-dark” gas in Perseus (Sec¬ 
tion ED, and finally summarize our conclusions (Section 

ED. 

2. BACKGROUND: METHODS TO ESTIMATE THE 
CORRECTION FOR HIGH OPTICAL DEPTH 

In most radio observations, HI is detected in emission, 
and the intensity of radiation is measured as the bright¬ 
ness temperature as a function of radial velocity, i.e., 
Tb{v). Since the HI optical depth (r) can be measured 
only via absorption line measurements in the direction of 
background radio continuum sources, the optically thin 
approximation of r 1 is frequently employed to esti¬ 
mate the HI column density: 

7V(HI) (cm-2) = 1.823 x j TB{v)dv (K km s"^). 

( 1 ) 
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Over the past three decades, several approaches have 
been employed to estimate how much of the true total 
HI column density is underestimated in the optically thin 
approximation. Most of these approaches can be classi¬ 
fied as “isothermal”. and th e only multiphase approaches 
are by iDickev et al.i (j200r)( l and HT03a. Here we sum¬ 
marize main results from some of the most important 
st udies. _ 

IDickev fc Benson! (I1982D used 47 emission/absorption 
spectral line pairs in the direction of background sources, 
and estimated the ratio of the HI column density from 
the absorption spectra to the HI column density in the 
optically thin approximation. Although ^1 at high and 
intermediate Galactic latitudes, the ratio reached ^1.8 at 
low latitudes. There was considerable scatter in the ra¬ 
tio, however: several lines of sight at low latitudes showed 
small ratios, suggesting that the low latitude directions 
with large ratios likely intersect dense molecular clouds. 
In order to compute the ratios, HI in each velocity chan¬ 
nel was assumed to have a single temperature (“isother¬ 
mal” ap proximation). 

In the IDickev et al.l (l2000l) study of the Small Mag¬ 
ellanic Cloud (SMC), HI absorption observations were 
obtained in the direction of 13 background radio contin¬ 
uum sources. The corresponding em ission spectra were 
derive d by averaging HI profiles from IStanimirovic et al.l 
()1999[1 over a 3 x 3 pixel region (pixel size = 30") centered 
on the position of each source. The correction factor for 
high optical depth was calculated for each velocity chan¬ 
nel in the isothermal approximation, and the line of sight 
integrated value / was expressed as a function of the un¬ 
corrected A^(HI): f = 1 + 0.667(logiQ A^(HI) — 21.4) for 
iV(HI) > cm“^. This relation was then applied to 

the 7V(HI) image of the SMC on a pixel-by-pixel basis, 
resulting in a ^10% increase of the total HI mass from 
^3.8 X 10® Mq to ~4.2 X 10® Mq. Although negligible 
at A^(HI) < 3x 10^^ cm“^, the correction factor increased 
with the uncorrected A^(HI) up to ^1.4 at iV(HI) ~ 10^^ 
cm“^. In some cases, the correction factors for individual 
channels were larger than the integrated value, reaching 
up to ~2. However, such values covered only a narrow 
range of channels, and their effect on A(HI) was rela¬ 
tively small. Finally, the authors rederived the correc¬ 
tion factor in the two-phase approximation, and found 
that the difference between the one- and two-phase cases 
depends on the relative location of cold and warm HI 
components along a line of sight ((Section l4^ for details). 

In the Millennium Arecibo 21-cm Absorption Line Sur¬ 
vey, HT03a obtained HI emission and absorption spectra 
toward 79 randomly positioned radio continuum sources, 
and performed Gaussian decomposition to estimate the 
physical properties of individual CNM and WNM com¬ 
ponents (column density, optical depth, spin tempera¬ 
ture Ts, etc.). These multiphase analyses showed that 
two or more components with very different spin tem¬ 
peratures can contribute to a single velocity channel, 
implying t hat the isotherma l treatm ent may not be sat¬ 
isfactory. iHeiles fc Trolandl (|2003bl) (HT03b hereafter) 
then calculated the correction factor using the Gaussian 
decomposition results, which they called i?raw = !//■ 
There were interesting variations in i?raw, ranging from 
~0.3 to ^1.0 (/ = ~1.0-3.0; Appendix B for details). 
Specifically, / ~ 1.3 was found for the Taurus/Perseus 
region. 


A v ery different approach was adopted in iBraun et al.l 
(|2009tl to calculate the correction for high optical depth 
in M31. They only used high-resolution HI emission 
observations for their modeling, and assumed that a 
single cold component determines the brightness tem¬ 
perature along a line of sight. While previous sim¬ 
ilar studies have applied a sin gle temperature to the 
images of entire ga laxies (e.g., iHenderson et al.l 1982; 
iBraun fc Walterb^ 1982). lSaun et alT i 200911 ^imated 
the spin temperature and non-thermal velocity disper¬ 
sion for each pixel. After excluding HI spectra that likely 
suffer from a high blending of different components along 
a line of sight, they noticed that the opaque HI is or¬ 
ganized into filamentary complexes and isolated clouds 
down to their resolution limit of ^100 pc. The spin tem¬ 
perature was found to increase from ^20 K to ^60 K with 
radius to 12 kpc, and then to decline smoothly down to 
~20 K beyond 25 kpc. The estimated correction resulted 
in a ^30% increase of the global HI ma ss of M31. Us¬ 
ing the same methodology, iBraunI (j2012tl found that the 
correction for high optical depth increases the HI masses 
of the Large Magellanic Cloud (LMC) and M33 by the 
same amount (^30 %). While the main advantage of the 
IBraun et al.l (|2009D approach is clearly that galactic-scale 
images of the opaque component can be produced solely 
from HI emission observations, the method has several 
weaknesses. For example, it does not consider multiple 
components along a line sight and how they self-absorb 
each other. In addition, it does not take account of the 
possibility that some of the brightness temperature could 
c ome from unabso r bing w arm HI components. 

IChengalur et al.l (|2013D tested the optically thin and 
isothermal approximations with Monte Carlo simulations 
of the multiphase ISM. They varied the fraction of gas in 
three phases (CNM, WNM, and the thermally unstable 
neutral medium) and the location of each phase along a 
line of sight. A wide range of values were assumed for 
the HI column density (10^°-10^'^ cm“^) and the spin 
temperature (20-5000 K). They found that the optically 
thin approximation underestimates the true HI column 
density by a factor of ^1.6 when J rdv ~ 1 km s“^, while 
the underestimate can be as high as a factor of ~20 when 
f rdv ^ 10 km s“^. On the other hand, the simulations 
showed that the isothermal estimate tracks the true HI 
column density to better than 10% even when f rdv ~ 5 
km s“^. Their conclusion that the isothermal estimate 
provides a good measure of the true HI column density of 
up to ^5 X 10^® cm“^ was insensitive to the assumed gas 
temperature distribution and the positions of the differ¬ 
ent p hases along a line of sig ht. We note that Equation 
(1) of IChengalur et ahl (|2ni3h does not include the con¬ 
tribution from the cosmic microwave background (CMB) 
and the Galactic synchrotron emission, which can be sig¬ 
nificant in certain cases (e.g., low Galactic latitudes). In 
addition, the authors did not consider self-absorption of 
the WN M by the foreground CNM. 

iLLsztl (j201^ compared A^( HI) from Galact i c HI surveys 
with E{B — V) derived by ISchlegel et al.l (|1998l) . and 
found a strong linear relation between fV(HI) and E{B — 
V) ^ 0.02-0.08 mag and a flattening of the relation at 
E{B — V)> 0.08 mag. While this flattening, likely due 
to H 2 formation , was essentially the same effect as what 
iLee et all (I2012D found for individual regions in Perseus, 
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Fig. 1. — 26 radio continuum sources overlaid on the HI column density image at 4' resolution (40+32.14 excluded; Section l3.ll for 
details). The HI column density image is produced by integrating the GALFA-HI cube from ulsr = to +15 km s~^, and the 

gray contours are from the CfA CO integrated intensity image at 8.4' resolution. The contour levels range from 10% to 90% of the peak 
value (69 K km s“^) with 10% steps. In this figure, both the Perseus and Taurus molecular clouds are seen. 


the relation derived bv iLisztl (l2014ll covers a large spatial 
area with randomly selected lines of sight predominantly 
at | 6 | > 20 ° By using HI absorption data compiled by 
iLiszt et "ni (|2010ll , the author then derived the correction 
for high optical depth, and applied it to the fV(HI) data. 
The flattening at E(B — V) > 0.08 mag persisted after 
the correction, confirming its origin in the onset of H 2 
formation. The derived correction factor increased from 
^1.0 at E{B — V) ^ 0.01 mag to ^1.4 at E{B — V) ^ 1 
mag, and was < 1.2 at E(B — V) < 0.5 mag. 

Recently, iFukui et'aP (|2015[1 suggested a new approach 
to estimate the correction for high optical depth by us¬ 
ing Planck dust continuum data. They noticed that the 
dust optical depth at 353 GHz (T 353 ) correlates with the 
HI column density, and the dispersion in this relation be¬ 
comes much smaller when the data points are segregated 
based on the dust temperature (Tdust)- The highest dust 
temperature was assumed to be associated with the opti¬ 
cally thin HI, and the saturation seen in the r 353 -V(HI) 
relation was then solely attributed to the optically thick 
HI. By coupling t he T 3 .t; 3 -jV(HI) rela tion with radiative 
transfer equations, iFukui etT] (120151) calculated Tg and 
T for the Galactic sky at | 6 | > 15° on a pixel-by-pixel ba¬ 
sis. They found that more than 70% of the data points 
have Tg < 40 K and r > 0.5, and similar results were 
obtained for the high latitude molecular clouds MBM 
53, 54, 55, and HLGG 92-35 (iFnkni et a.l.l[2(rn . The 
correction for high optical depth resulted in a factor of 
^2 increase in the total HI mass in the solar neighbor¬ 
hood, implying that the optically thick HI may explain 
the “GO-dark” gas in the Galaxy. 

3. DATA 

3.1. HI Emission and Absorption Observations 

We use the HI emission and absorption observations 
from Paper I. The observations were performed with 
the Arecibo telescope^ using the L-band wide receiver, 
and were made toward 27 radio continuum sources lo- 

^ The Arecibo Observatory is operated by SRI International un¬ 
der a cooperative agreement with the National Science Foundation 
(AST-1100968), and in alliance with Ana G. Mendez-Universidad 
Metropolitana and the Universities Space Research Association. 


cated behind Perseus. The target sources were se lected 
from the NRAO VLA Sky Survey (|Condon et al.l 1998) 
based on flux densities at 1.4 GHz greater than ^0.8 
Jy, and are distributed over a large area of ~500 deg^ 
centered on the cloud (Figure [1])®. The angular resolu¬ 
tion of the Arecibo telescope at 1.4 GHz is 3.5'. For 
the observations, a special procedure was adopted to 
make a “17-point pattern”, which includes 1 on-source 
measurement and 16 off-source measurements (HT03a; 
IStanimirovic fc HeilesI 2005). This procedure was de¬ 
signed to consider HI intensity variations across the sky 
and instrumental effects involving telescope gains. The 
data were processed using the reduction software devel¬ 
oped by HT03a, and the final products for each source in¬ 
clude an HI absorption spectrum an “expected” 

HI emission spectrum (Texp(n); HI profile that we would 
observe at the source position if the continuum source 
were not present), and their uncertainty profiles. Among 
the 27 sources, 4G-I-32.14 was excluded from further anal¬ 
yses because of its saturated absorption spectrum. With 
an average integration time of 1 hour, the root-mean- 
square (rms) noise level in the optical depth profiles was 
^1 X 10“® per 1 km velocity channel. Finally, the de¬ 
rived optical depth and “expected” emission spectra were 
decomposed into separate CNM and WNM components 
using the technique of HT03a, and physical properties 
(optical depth, spin temperature, column density, etc.) 
were computed for the individual components. We refer 
to Paper I for details on the observations, data reduction, 
line fitting®, and GNM/WNM properties. 

3.2. HI Emission Data from the GALFA-HI Survey 

In order to evaluate different methods for deriv¬ 
ing the correction for high optical depth, we also 

® In this paper, we quote all velocities in the local standard of 
rest (LSR) frame, which is defined based on the average velocity 
of stars in the solar neighborhood: 20 km s“’ toward (R.A.,deck) 
= (18*',30°) in B1900. 

® Pros and cons of our Gaussian fitting method were discussed 
in HT03a in detail. In the future, we plan to compare the Gaussian 
fitting method with results from numerical simulations to investi¬ 
gate biases that could be introduced by Gaussian fitting (Lindner 
et al. in prep). 
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use the HI emission dat a from the GALFA-HI survey 
(|Stanimirovic et ^ 2006: [Peek et al.l 2011). GALFA-Hl 
uses ALFA, a seven-beam array of receivers at the focal 
plane of the Arecibo telescope, to map the HI emission in 
the Galaxy. Each of the seven dual polarization beams 
has an effective beam size of 3.9 ' x 4.1'. 

For Perseus, iLee et al] (|2012il produced an HI cube 
centered at (R.A.,decl.) = (03‘'29“52",-b30°34'l") in 
J2000^° with a size of ~15° x 9° by combining a num¬ 
ber of individual GALFA-Hl projects. We use the same 
data here, but extend the HI cube up to ^60° x 18° to 
include all radio continuum sources in Paper 1. The HI 
column density image derived from the extended HI cube 
is shown in Figure 1 along with our continuum sources 
(4C-I-32.14 excluded). 



3.3. HI and H 2 Distributions of Perseus 


We use the fV(Hl) and 7V(H2) images from iLee et al.l 
(|2012ll . To derive the iV(HI) image, we integrated the HI 
emission from ulsr = km s“^ to -1-15 km s“^ in the 
optically thin approximation. This velocity range was de¬ 
termined based on the maximum correlation between the 
A^(Hl) image and 2MASS-based Ay data from the GOM- 
PLETE Survey (iRidge et al.ll2006ll . In order to construct 
the N{R 2 ) image of Perseus with a large sky coverage, 
we use d the 60 /zm and 100 iim da ta from the IRIS 
survey (|Miville-Deschenes fc Lagachel 20051. and derived 
the dust optical depth at 100 /rm (rioo) by assuming that 
dust grains are in thermal equilibrium. For this purpose, 
the emissivity spectral index of ,5 = 2 was adopted, and 
the contribution from very small grains (VSGs) to the 
intensity at 60 fim (/eo) was removed by calibrating the 
derived TH„st image with DlRBE-based Tdust data from 
ISchlegel et al.l (jl998ll . We then converted the tiqo image 
into the Ay image by finding the conversion factor X 
for Ay = Atioo that results in the minimum difference 
between the derived Ay and the COMPLETE Ay. This 
cal ibration of rmn t o the C OMPLETE Ay was motivated 
bv [Goodman et ahi (l2009ll who showed that dust extinc¬ 
tion at near-infrared (NIR) wavelengths is the best probe 
of total gas column density. Finally, we measured a local 
dust-to-gas ratio (D/G) by examining the Av-A^(Hl) re¬ 
lation for diffuse regions, and derived the A^(H 2 ) image 

by 


7V(H2) = - 


Ay 

D/G 


- N{m) 


( 2 ) 


The derived A^(Hl) and iV(H 2 ) images are at 4.3' resolu¬ 
tion (corresponding to ^0.4 pc at the distance of 300 pc), 
and their median Icr uncertainties are ^5.6 x 10^^ cm~^ 
and_~3.6 x 10^® cm“^. See Sections 3 and 4 of ILee et'al] 
(|2012ll for details on the derivation of the A^(Hl) and 
A^(H 2 ) images and their uncertainties. 


3.4. CO Data from the CfA Survey 

We use CO in t egrated intensity (Ian) data from 
iDame et al.l (l200lD . iDame et al.l (1200111 produced a com¬ 
posite CO survey of the Galaxy at 8.4' resolution by com¬ 
bining individual observations of the Galactic plane and 
local molecular clouds. The observations were conducted 
with the Harvard-Smithsonian Center for Astrophysics 


Fig. 2.— Comparison of the HI column densities calculated using 
the two different methods: A^exp from the spatial derivative method 
versus A^cube from the simple averaging method. Both quantities 
are estimated in the optically thin approximation, and the black 
solid line shows a or^±o-o^e relation. ^ 

(CiA) telescope, ihe rinai cube had a unirorm rms noise 
of 0.25 K per 0.6 5 km s~^ vel o city c hannel. To estimate 
Ico for Perseus, iDame et al.l (l2001ll integrated the CO 
emission fro m ut,sr = —15 km s“^ to -1-15 km s“^. See 
Section 2 of iDame et al.l (|2001ll for details on the obser¬ 
vations, data reduction, and analyses. 

4. CORRECTION FOR HIGH OPTICAL DEPTH 

In the direction of 26 radio continuum sources, we mea¬ 
sured the optical depth profiles which we use to estimate 
the true total HI column density, iVtot- Along the same 
lines of sight, we also have the emission spectra that can 
be used to calculate the HI column density in the opti¬ 
cally thin approximation, Viow-r- This column density 
would be the only available information if no HI absorp¬ 
tion data were present. In this section, we examine how 
/ = IVtot/Viow-T, which we call the correction factor for 
high optical depth, varies with Aiow-r- Our aims are to 
offer an analytic estimate of /(Mow-t) for Perseus us¬ 
ing our HI emission and absorption measurements, and 
then to a p ply th is correction to the V(HI) image from 
ILee et aP (|20I2ll on a pixel-by-pixel basis. In this way, 
we can account for the optically thick HI that was missed 
in the GALFA-HI emission observations. Since there are 
several approaches to derive /(Aiow-t), we first compare 
different methods using full line of sight information. 

4.1. Calculating the HI Column Density in the 
Optically Thin Approximation (N\ov,-t) 

As the HI emission spectrum in the direction of a radio 
continuum source is affected by absorption, it is not pos¬ 
sible to obtain an emission profile along exactly the same 
line of sight as probed by the HI absorption observation. 
For this reason, we instead estimated the “expected” HI 
emission spectrum Texp(u), which is the profile we would 
observe if the continuum source turned off, by modeling 
the “17-point pattern” measurements. In this modeling, 
spatial derivatives of the HI emission (up to the second 
order) were carefully taken into account (Section 2.1 of 
Paper I for details). Then the HI column density in the 
optically thin approximation, N^xp, can be calculated by 

iVexp (cm“^) = 1.823 X J Texp{v)dv (K km s“^). 

( 3 ) 


In this paper, we quote all coordinates in J2000. 
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Fig. 3. — METHOD 1: (left) / = A^tot/^exp (A^tot: derived using the Gaussian decomposition results) as a function of logio(Afexp/10^^). 
The blue squares show the (l/(7^)-weighted mean values in 0.2-wide bins in logio(Afexp/10^^), and the linear fit determined for all 26 data 
points is indicated as the green solid line (Equation [6)1 . (right) / as a function of the integrated optical depth. The (l/(T^)-weighted mean 
values in 3.4 km s“^-wide bins in the integrated optical depth are presented as the blue squares, and the green solid line shows the linear 
fit to all 26 data points: / = (0.024 ± 0.004) J r{v)dv + (1.019 ± 0.019). 


We compute A^exp over the velocity range where Texp(u) 
is above its 3cr noise level. This velocity range covers all 
CNM and WNM components for each source, and the 
median velocity range for all 26 sources is ulsr = “39 km 
s“^ to -1-20 km s“^. We then estimate the uncertainty in 
Nexp by propagating the Texp(ii) error spectrum through 
Equation ([3]), finding a median of ~3.0 x 10^® cm“^. 

Additionally, we use spectra from the GALFA-HI cube 
(pixel size = 1') in the direction of our radio continuum 
sources to derive the “expected” emission profiles. This 
simpler approach has been employed when HI absorp¬ 
tion spectra were obtained without any special str ategy 
such as the “17-point pa ttern” fe.g..lDickey et al.l 2000; 
iMcClure-Griffiths et al.l 2001: IPickev et al.l 2003). We 
extract HI emission spectra from a 9 x 9 pixel region 
(roughly the “17-point pattern” grid size; Figure 1 of 
HT03a) centered on each continuum source. As the HI 
spectra right around the continuum source are likely af¬ 
fected by absorption of the background emission, we ex¬ 
clude the HI spectra from the central 3x3 pixel region 
(roughly one Arecibo beamwidth across). By averaging 
the remaining 72 spectra, we then compute an average 
emission spectrum Tavg(u), and estimate the correspond¬ 
ing HI column density fVcube by 

A^cube (cm“^) = 1.823 x J Ts,vg{v)dv (K km s“^). 

. ( 4 ) 

The uncertainty in A^cube is estimated by calculating the 
standard deviation of the extracted 72 spectra, crT&vgiv), 
and propagating it through Equation ([4]). The median 
value for all 26 sources is ~1.4 x 10^® cm“^. 

In Figure [2l we compare the HI column densities cal¬ 
culated using the two methods. Most sources probe the 
HI column density of ~(5-16) x 10^° cm“^, and the 
last five sources (3C132, 3C131, 3C133, 4C-I-27.14, and 
4C-I-33.10) extend this range by a factor of ^2. We find 
an excellent agreement between the two methods up to 
^3 X 10^^ cm“^. The ratio of A^exp to A"cube ranges from 
~0.9 to ~1.1 with a median of ^1.0, suggesting that the 
two HI column density estimates are consistent within 
10%. Considering the more careful examination of spa¬ 
tial variations in the HI emission, we continue by using 


Nexp as A^iow-r in the following sections. 

4.2. METHOD 1 - Gaussian Decomposition to 
Estimate A^tot 

In Paper I, we performed Gaussian decomposition of 
the optical depth and “expected” emission spectra, and 
calculated the properties of individual CNM and WNM 
components while considering self-absorption of both the 
CNM and the WNM by the CNM. All Gaussian decom¬ 
position results (peak brightness temperature, peak op¬ 
tical depth, spin temperature, etc.) for each component 
are presented in Table 2 of Paper I. These results enable 
us to derive the true total HI column density along a line 
of sight by 

Wot (cm“^) = A^cnm + Awnm 

= 1.823 X 10^® f (^ 

0 

K-1 

-f (K km s"^), 

0 

(5) 

where the components with subscript n refer to the 
CNM, the components with subscript k refer to the 
WNM, To is the peak optical depth, vq is the central 
velocity, Tq is the peak brightness temperature, and Sv 
is the 1/e width of the component. Here Wot is cal¬ 
culated over the velocity range determined as having 
Fexp(^^) higher than its 3cr noise. For the uncertainty in 
Wotj we use errors of the fitted parameters provided by 
Gaussian decomposition to perform a Monte Carlo simu¬ 
lation where 1000 A"cnm and A^wnm values are computed 
from normally distributed parameters. The standard de¬ 
viations of the A'cnm and A^wnm distributions are then 
added in quadrature to estimate the uncertainty in Wot • 
This method of deriving Wot was used by HT03b for 
their HI absorption measurements toward background 
sources randomly located over the whole Arecibo sky. In 
this study, we focus on a localized group of background 
sources in the direction of Perseus. 
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Fig. 4.— METHOD 2: (left) Correction factor per velocity channel /chan('^) (estimated in the isothermal approximation) for all 26 
sources, (right) Histogram of /chan('^) values. 


The (integrated) correction factor, / = Vtot/Mow-r = 
-Vtot/Vexp, is shown in Figure |3] (left) as a function of 
iVexp- Clearly, the correction factor increases with Ng-^p- 
We then present the (l/CT^)-weighted mean values as the 
blue squares and the linear fit to all 26 data points as 
the green soild line^^: 

/ = logio(Vexp/10^°) xa + b 

= logio(Vexp/102°)(0.32 ± 0.06) + (0.81 ± 0.05). 

In general, the correction factor ranges from ~1.0 at ~3.9 
X 10^° cm“^ to ~1.2 at ^1.3 x 10^^ cm“^ (maximum 
uncorrected HI column density in Perseus). While / and 
fVexp show a good correlation (Spearman’s rank correla¬ 
tion coefficient of 0.80), there are two sources with rela¬ 
tively high correction factors at ^10^^ cm“^, 3C092 and 
3C093.1. Interestingly, the two are located behind the 
main body of Perseus (Figure 1). Their high / values of 
^1.5-1.6 could result from an increased amount of the 
cold HI in the molecular cloud relative to the surround¬ 
ing diffuse ISM. The CNM fraction is indeed ^0.4 for 
both sources, which is higher than the median value of 
~0.3 for all 26 sources (Section 4.3 of Paper I). However, 
this is not the maximum CNM fraction in our measure¬ 
ments (~0.6). Observing a denser grid of radio contin¬ 
uum sources behind Perseus and repeating the calcula¬ 
tions would be an interesting way to test the cold HI hy¬ 
pothesis. Finally, the correction factor is also presented 
as a function of the integrated optical depth in Figure |3] 
(right). As expected, there is a clear correlation (Spear¬ 
man’s rank correlation coefficient of 0.94). 

We note that our results are not sensitive to HI com¬ 
ponents at wlsr < ~20 km s“^, which are likely unasso¬ 
ciated with Perseus (Section 4.2 of Paper I); limiting the 


calculation of Vexp and Wot to wlsr > — 20 km s“^ or 
excluding the five sources showing such HI components 
at large negative velocities (corresponding to the sources 
with logio(Wxp/10^°) > 1.2 in Figured 3C132, 3C131, 
3C133, 4C-I-27.14, and 4C-I-33.10) results in linear fit co¬ 
efficients that are consistent with what we present here 
within uncertainties.^^ 


4.3. METHOD 2 - Isothermal Estimate o/Wot 


By assuming that eac h velocity channel r epresents gas 
at a single temperature, iDickev et al.l (|2000ll showed that 
the correction factor per velocity channel can be written 
as 


/chan('^^) 


CoTs{v)t{v) 
CoTexp(n) ’ 


(7) 


where Co = 1.823 x 10^® cm ^/(K km s ^). In addition, 
Texp{v) was expressed as 

Tg^p{v) = Tg{v){l-e-^^^^). ( 8 ) 

This equation assumes the absence of any radio contin¬ 
uum source behind the absorbing HI cloud. As a result, 
fcha.niv) simply becomes 


fchani’^'j 


t{v) 

1 - ■ 


(9) 


However, while the radio continuum source is absent in 
Equation ([S]), some diffuse radio continuum emission is 
always present, and should not be ignored. This emission 
includes the CMB and the Galactic synchrotron emission 
that varies across the sky and becomes strong toward the 
Galactic plane. We call a combination of these contribu¬ 
tions as Tsky, and Equation ([5]) then has to be rewritten 
as 


In attempting to be consistent with our calcul ation of the un¬ 
certainty in / for the isothermal method ^Section I4.3II . we run a 
full Monte Carlo simulation based on errors of the fitted parame¬ 
ters from Gaussian decomposition. In this simulation, 1000 Nexp 
and A^tot values are calculated from normally distributed param¬ 
eters, and the standard deviation of 1000 / values is used as the 
uncertainty in /. We find that linear fit results from using this 
error estimate, / = logio(A^ex^lO^*^)(0.25 ± 0.03) -h (0.87 ± 0.02), 
are consistent with Equation {Sjl within uncertainties. 


TLpiv) = Tgiv)il - e-"('’)) + Tskye-^^’’). (10) 

To be specific, limiting the calculation of A^exp and A^tot to 
'I’LSR > —20 km results in / = logio(Afexp/10^°)(0.44±0.07)-|- 
(0.73 ± 0.06) and / = (0.034 ± 0.005) / T(v)dv + (1.004 ± 0.020). 
Similarly, excluding the five sources that show the HI components 
at ulsr < —20 km s”^ leads to / = logio(Alexp/10'^'^)(0.31±0.11)-|- 
(0.82 ± 0.09) and / = (0.040 ± 0.011) / T{v)dv + (0.988 ± 0.027). 
































Fig. 5.— METHOD 2: (left) / = A^tot/-^exp (-^tot- estimated in the isothermal approximation) as a function of logio(A^exp/10^^). The 
blue squares show the (l/(T^)-weighted mean values in 0.2-wide bins in logio(A^exp/10^®), and the linear fit determined for all 26 data 
points is indicated as the green solid line (EQuation ll5|l . (right) / as a function of the integrated optical depth. The (l/cr^)-weighted mean 
values in 3.4 km s“^-wide bins in the integrated optical depth are presented as the blue squares, and the green solid line shows the linear 
fit to all 26 data points: / = (0.023 zb 0.001) J T(v)dv -|- (1.030 ib 0.006). 


Considering that HI emission spectra are generally base¬ 
line subtracted during the reduction process, Tgky can be 
removed from both sides of Equation (Unj. Thenrexp(^;) 
= ’^e*xp('*^) —"^sky, the quantity we have been working with 
so far, can be expressed as 

Texpiv) = r'exp(l^) ~ Tsky 

= (T,(n)-r,ky)(l-e-"(’')). 

As a consequence, the correction factor becomes 
r ^ ^ _ Tsiv) t{v) 

/chanW 7;(r;)-rsky 1-e-W’ 

or with direct observables, 

r ^ ^ '^(^) , 

/chanW-Jskyy^^^^^^ ^ 1 - ' 

In order to estimate the contribu tion from the Galac - 
tic synchrotron emission, we use the lHaslam et al.l (|1982[1 
408 MHz survey of the Galaxy. The brightness tempera¬ 
ture at 408 MHz is converted to 1.4 GHz using the spec¬ 
tral index of —2.7. As the absolute Galactic latitude of 
our continuum sources is generally higher than 10°, the 
synchrotron contribution is small with T^ky ranging from 
2.78 K to 2.80 K (Table 1 of Paper I). Based on the his¬ 
togram of Ts for the individual CNM components (Figure 
5b of Paper I), we can then provide a rough estimate of 
Ts{v)/[Ts{v) — Tsky]: the expected range is narrow, from 
^1.0 to ~1.2. Clearly, for molecular clouds located closer 
to the Galactic plane the contribution from the diffuse 
radio continuum emission will be more significant. 

In Equation ([7]), /chan(^^) essentially represents the cor¬ 
rection that needs to be applied to Texp(r’) to calculate 
the true brightness temperature profile. As a result, the 
true total HI column density can be obtained by 


( 11 ) 

( 12 ) 

(13) 


Vtot (cm ^) = 1.823x10^® / fchs.niv)Texpiv)dv(Kkms ^). 

(14) 

We derive /chan(^^) for all 26 sources (Equation (TS]), and 
present the results in EigurelH About 91% of the /chan('c) 
values are between 1 and 2, and the fraction of velocity 


channels with /chan(i^) > 2 is very small. We then calcu¬ 
late Wot using Equation da, and show the (integrated) 
correction factor, / = Wot/Wxp; as a function of N^^p 
in EigureO (left). Here the integration is done over the 
velocity range where T'exp(n) is higher than its 3cr noise. 
The uncertainty in / is derived by running a Monte Garlo 
simulation where the optical depth and “expected” emis¬ 
sion error spectra are propagated through Equations m 
and m to compute 1000 / values. The standard devia¬ 
tion of the / distribution is used as the final uncertainty 
in /. 

Similar to the Gaussian decomposition method, we find 
a good correlation between / and Wxp (Spearman’s rank 
correlation coefficient of 0.84). The linear fit determined 
using all 26 data points is 

/ = logio(Wxp/102°)(0.25 ± 0.02) + (0.87 ± 0.02). (15) 

Additionally, / is plotted as a function of the integrated 
optical depth in Eigure[5] (right), again showing a clear 
correlation (Spearman’s rank correlation coefficient of 
0.97). 

Both graphs in Eigure[5]are very similar with those in 
Figure[3]for the Gaussian decomposition method. Specif¬ 
ically, the linear fit coefficients are consistent within un¬ 
certainties. This is surprising considering that the two 
methods are very different. In particular, the isothermal 
method assigns a single spin temperature to each velocity 
channel, while the Gaussian decomposition method al¬ 
lows a single velocity channel to have contributions from 
sever al HI compon e nts wit h different spin temperatures. 

In iDickev et al.l (j200r)( l. the authors updated the 
isothermal method by incorporating the two-phase ap¬ 
proximation. As input parameters, this method then 
required the spin temperature of the cold HI and the 
fraction of the warm HI that is in front of the c o ld HI , 
the quantity they referred to as q. iDickev et al.l (l200Clf l 
showed that for their SMG data there is no difference be¬ 
tween the one- and two-phase approximations regarding 
the correction factor if q > 0.5, while the difference be¬ 
comes more pronounced when q < 0.25. In the Gaussian 
decomposition method, this fraction (Fk in Paper I) is 
important as well, but is difficult to constrain. Thus, the 
htting process was repeated for Fk = 0, 0.5, and 1, and 
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Correction Factor 
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Fig. 6. — Corr ection factor at 4.3' resolution estimated in Section l4.2l (Equation I51l. In addition, the 3cr contour of the “old” fV(H 2 ) from 
ILee et aP II2012I 1 (before the correction for high optical depth) is overlaid in gray, while the 3cr contour of the CfA Iqq is shown in black. 
The resolutions of the ^(112) and Iqq images are 4.3' and 8.4'. 


these results were used to estimate the final uncertainty 
in Ts (Section 3 of Paper I and HT03a for details). 

To understand why the two different methods result 
in similar correction factors, we compare Equations m 
and (fT4|) (Appendix A), and find that they become com¬ 
parable regardless of Fk when r <C 1 and Tgky Tg. 
In our observations of Perseus, the median peak opti¬ 
cal depth for all CNM components is ^0.2, and only 
a small number of the components has the peak opti¬ 
cal depth higher than 1 (10 out of 107; Section 4.1 of 
Paper I), satisfying the condition. In addition, we al¬ 
ready showed that Tgky is small with ~2.8 K for Perseus. 
The difference between the Gaussian decomposition and 
isothermal methods, however, will be more significant for 
molecular clouds that have a large amount of the cold, 
optically thick HI and/or a substantial contribution from 
the diffuse radio continuum emission. Due to the more 
self-consistent way to derive Texp(u), we continue by us¬ 
ing the Gaussian decomposition results for further anal¬ 
yses. 

Finally, we note that Equations (O and m could be 
biased against very high optical depths as they result in 
saturated absorption spectra, e.g., 4G-I-32.14, the source 
we had to exclude from our analyses due to its highly sat¬ 
urated absorption profile. In addition, the equations are 
based on our explicit assumption of the linear relation be¬ 
tween / and logio(lVexp/10^°). The HI emission and ab¬ 
sorption measurements obtained by HT03a,b along many 
random lines of sight through the Galaxy show that the 
linear relation indeed describes the observations well up 
to logio(fVexp/10^°) ~ 1.5, and the fitted coefficients are 
consistent with Equations ([6]) and (1151) within uncertain¬ 
ties (Appendix B). There is some interesting deviation 
from the linear relation, however, particularly for a few 
sources with logio(Aiexp/10^°) > 1. This deviation could 
suggest a non-linear relation at high column densities, 
and its significance needs to be further examined with 


more HI absorption measurements. In the future, it will 
also be important to study the dense CNM using alterna- 
tive tracers s uch as Cl and CH fine-structure lines (e.g., 
iPineda et al.l 2013). 

4.4. Comparison with Previous Studies 

Dickey et al. (2000 ) . Th e correction factor calcu¬ 
lated by iDickev et al.l ()2000D for the SMC using full 
line of sight information can be rewritten as / = 
logio(iVexp/102°)0.667 -b 0.066 for iVexp > 2.5 x lO^i 
cm“^. The range of the HI column density we probe fo r 
Perseus barely overlaps with that in iDickev et al.l (|2000D . 
as the HI column density in the low-metallicity SMC 
is significantly hig her compared to the Galaxy. While 
the difference with iDickev et al.l (|2000D in / values de¬ 
pends on i Vpxni our correction f actor is only ^4% higher 
than what IDickev et al.l (j200r)D suggests when extrapo¬ 
lated to t heir maximum HI co lumn density of 10^^ cm“^. 
Similarly, IDickev et al.l (120031) used HI emission and ab- 
sorp tion data from the Sou thern Galactic Plane Sur¬ 
vey (|McChire-Griffiths et al.l 2005) in combination with 
the isothermal method, and estimated / = ~1.4-1.6 for 
sources located at 326° < I < 333° and |6| < 1°. This 
correction factor is comparable with what we find for 
Perseus. 

Heiles & Troland (2003a,b). HT03a,b performed Gaus¬ 
sian decomposition of 79 HI emission/absorption spectral 
line pairs, and derived the correction factor, which they 
called i?raw = !//• They found that i?raw ranges from 
^0.3 to ~1.0, corresponding to / = ^1.0-3.0 (Appendix 
B). In particular, they estimated / ^ 1.3 for the Tau¬ 
rus/Perseus region, which is similar with what we find 
for Perseus. 

Braun e t al. (2009). Our corre ction factor i s smaller 
than what iBraun et al.l (j2009t) and lBraiii] (j2012t) claimed 
for M31, M33, and the LMC based on the modeling of 
HI emission spectra. They found that the correction ex- 
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Fig. 7.— Corrected ^(111) image at 4.3' resolution. 


ceeds an order of magnitude in many cases, and increases 
the global HI mass by ~30%. Even when considering the 
correction factor per velocity channel, we find the maxi¬ 
mum /chan(i^) 4 for only one source and /chan(u) < 3 
for the rest of our sources. 

Chengalur et al. (2013). Our corr ection factor can 
also b e compared with predictions from IChengalur et al.l 
(1201,31 1. This study performed Monte Carlo simulations 
where observationally motivated input parameters such 
as the column density, the spin temperatures of the CNM 
and the WNM, and the fraction of gas in each of the 
different phases were provided for ISM models. Our cor¬ 
rection factor versus integrated optical depth plots. Fig¬ 
ures [3] (right) an d [5] (right), can be dire ctly compared 
with Figure lA in lC^engalur et all (1201,^. We find that 
the correction factor bv IChengalur et al.l (|2013[1 is signif¬ 
icantly higher than our estimate, although the general 
trend of increasing correction wi th the integrated optica l 
depth is similar. For example, IChengalur et al.l (l201.3f l 
expects / ^ 20 when f rdv ~ 10 km s“^, while we find 
only / < 1.5. Similarly, our Figures [3l ( left) and [5] (left) 
can b e compared with Figure 2A in IChengalur et al.l 
(|2013ll . We find tha t our estimate is c onsist ent with the 
correction factor by IChengalur et al.l (1201311 for the col¬ 
umn d ensity less than 10^^ cm“^, while IChengalur et al.l 
(j201.3ll overestimates at the high end of our column den¬ 
sity range. If we extrapolate our relation to 10^^ cm“^, 
we expect a ^10 t imes lower correction factor than what 
IChengalur et al.l (j201,3ll suggests. The reason for their 
very high correction factor could be the inclusion of ex¬ 
tremely high column densities (lO^^-lO^'^ cm“^) in their 
ISM models, although it is not clear why it less affects 
the isothermal estimate of the HI column density. In the 
simulations, the median ratio of the true HI column den¬ 
sity to the isothermal estimate was ~1.1 when f rdv ~ 
5 km s“^. 

Liszt (2014) . Usin g the HI abs orption data compiled 
by iLiszt et ^ (|2010ll , iLiszd (l2014ll estimated the correc- 



Fig. 8.— (left) Normalized histograms of the two Al(HI) images, 
before (black) and after (gray) the correction for high optical depth. 
The median of each histogram is shown as the dashed line, (right) 
Same as the left panel, but for Af(H 2 ). 

tion factor for radio continuum sources located at high 
Galactic latitudes. While they did not provide detailed 
information about how exactly the derivation was done, 
their correction factor was small with / less than 1.2 for 
E{B — F) < 0.5 mag. This is compa rable to our f i nding . 

Fukui et al. (2014,2015). Finally. iFukui et al.l (|20I5 ) 
estimated the correction factor for the Galaxy at |5| > 15° 
by exploring the relation between rasa and N (HI) at 33' 
resolution. Their Figure 13 shows that the correction fac¬ 
tor distribution ranges from ^I.O to '^3.0, and p e aks at 
~1.5. Using the same methodology, iFukui etM] (|20I4D 
found that the correction increases the total HI mass of 
the high latitude molecular clouds MBM 53, 54, 55, and 
HLCG 92-35 by a factor of ^2. In general, the correction 
factor by Fukui et al. (2014,2015) appears higher than 
our estimate for Perseus. While we do not perform a de¬ 
tailed comparison with Fukui et al. (2014,2015), we test 
their claim of the optically thick HI as an alternative of 
the “CO-dark” gas in Section [T] 
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Fig. 9. — Rederived A^(H 2 ) image. The blank pixels correspond to regions with possible contaminations (point sources, the Taurus 
molecular cloud, and the “warm dust ring”; Section 4 of ILee et al.l 2012 for details). The 3cr contour of the new A^(H 2 ) is overlaid in 
gray, while the Scr contour of the CfA Iqq is shown in black. The resolutions of the N(H.2) and Iqq images are 4.3^ and 8.4'. The red 
rectangular boxes show the boundaries of the selected dark (B5, BIE, and Bl) and star-forming (IC348 and NGC1333) regions. 


5. APPLYING THE CORREGTION FOR HIGH OPTIGAL 
DEPTH TO PERSEUS 

We apply the correction de rived in S e ction 14.21 to the 
V(HI) image of Perseus from iLee et alJ (|2012f) for pixels 
with logio(-/Vexp/10^°) > 0.6 (where the correction factor 
is higher than 1). The correction factor and corrected 
V(HI) images are shown in Figures H] and 0 In addition, 
we present the normalized histograms of the two iV(HI) 
images, before and after the correction, in Figure El 
As Figure El indicates, the correction does not make a 
significant change in V(HI). To be specific, the median 
V(HI) increases by a factor of ~1.1 from ^7.9 x 10^° 
cm“^ to ~8.7 X 10^° cm“^, while the maximum V(HI) 
increases by a factor of ~1.2 from ^1.3 x 10^^ cm“^ 
to ~1.6 X 10^^ cm“^. In terms of the total HI mass, 
the correction results in a ^10% increase from ~2.3 x 
10^ Mq to ~2.5 X 10"^ Mf7i . This increase in the HI 
mass is comparable to what iDickev et al.l (|2000[1 found 
for the SMC, but is small er than the v alue estimated 
bv iBraun et al.l (1200911 and iBraiml (1201211 for M31, M33, 
and the LMC. In addition, we note that our correction 
factor image looks smooth compared to a gr anulated ap- 

S ice o f the co r rected V(HI) images bv iBraun et al.l 
and iBrauiil (1201211 . although Perseus (~80 pc x 
50 pc) would be unresolved or only marginally resolved 
in their studies. Finally, the HI mass increase due to the 
optically thick gas in Perseus is smaller than what Fukui 
et al. (2014,2015) derived for Galactic molecular clouds. 

6. HOW DOES THE HIGH OPTICAL DEPTH AFFECT THE 
HI SATURATION IN PERSEUS? 

6.1. Rederiving V(H 2 ) 

To investigate the impact of high optical depth on the 
observed HI saturation in Perseus, we first rederive the 
V(H 2 ) image using the correc ted V(HI) . In es sence, we 
use the same methodology as iLee et aO (|2012[1 : the Ay 


image is derived using the IRIS 60/100 fim and 2MASS 
Ay data, and a local D/ G is adopted to es timate V(H 2 ). 
We refer to Section 4 of iLee et al.l (I2012D for details on 
the method for deriving V(H 2 ) and its limitation, and 
summarize here main results. 

1. Gontamination from VSGs: To estimate Tjust from 
the ratio of leo to Iiooj the contribution from stochasti¬ 
cally heated VSGs to Iqq must be removed. For this pur¬ 
pose, we compar e our IRIS-based T^ns t with the DIRBE- 
based Tdust from ISchlegel et all (Il998ll . and find that the 
contrib ution from VSGs to Jeo is 78%. This is the same 
as what iLee et aH (I2012D found. 

2. Zero point calibration for tioq: We refine the zero 
point of the rioo image by assuming that the dust col¬ 
umn density traced by tiqo is proportional to V(HI) for 
atomic-dominated regions. Based on the zero point of the 
tioo“A^(III) relation, we add 1.1 x 10~ ^ to the r i nn im - 
age. This is slightly smaller than what ILee et aH (120121) 
added, i.e., 1.8 x 10“"^. 

3. Conversion from rioo to Ay: We convert rioo into 
Ay by adopting X — 740 for Ay = ATtioo that results 
in the best agreement betw een our IRIS-based Ay and 
the 2MASS-based Ay from lRidee et al.l (I2006D. This i s 
slightly higher than X = 720 used bv ILee et al.l (l2012ll . 
We compare our rederived Ay with the Ay image from 
ILee et aH (|2012D . and find that the ratio of the new Ay 
to the old Ay ranges from ^0.93 to ^1.02. 

4. Deriving a local D/G and V(H 2 ): We examine the 
Av-V(HI) relation, and find that the slope of Ay /V(HI) 
= 1 X 10“^^ mag cm^ is a good measure o f D/G for 
Perseu s. This is slightly lower than what ILee et aP 
(|2012f) estimated, i.e., 1.1 x 10“^^ mag cm^, and makes 
sense considering that our re derived Ay is esse ntially the 
same with the Ay image from iLee et al.l (|2012f) . while the 
corrected V(HI) is slightly higher than the uncorrected 
V(HI). Finally, we derive V(H 2 ) using Equation (E]), and 
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Fig. 10.— (left) Sjji versus Ejji + Sh 2 for IC348 and BIE (this study). All finite pixels in the rectangular boxes shown in Figure|9]are 
used for pl ottin g. The median Scr values of Shi and Shi + Sh 2 are indicated as the black dashed lines, while the best-fit curves d etermined 
in Section 16.41 are overlaid in red. The best-fit parameters are shown in the top right corner of each plot, (right) Plots from [Lee et al.l 


mask pixels with possible contaminations (point sources, 
the Taurus mo lecular cloud, and the “warm dust ring”), 
following what iLee et aP (|2012D did. We show the final 
Af(H 2 ) image in Figured as well as the normalized his- 
tograms of the new (this study) and old ()Lee et al.ll2012j) 
Af(H 2 ) images in Figure |8] (right). 

Conclusion: In our rederivation of Af(H 2 ), all param¬ 
eters are identical w ith or only slightly different from 
what iLee et al.l (|2012h used. As a result, the new Af(H 2 ) 
is comparable to the old N(ll 2 ), shown as a good agree¬ 
ment between the two histograms in Figure [8] (right). 
The rederived N{li 2 ) ranges from —1.5 x 10^° cm“^ to 
5.1 X 10^^ cm“^ with a median of 4.6 x 10^® cm“®, and 
~83% of the pixels whose S/N is higher than 1 have the 
new Af(H 2 ) differing from the old Af(H 2 ) by only 10%. 
Finally, we note that ^30% of all finite pixels have nega¬ 
tive Af(H 2 ) values of mostly around —(1-5) x 10^® cm“®, 
which are very close to zero considering the median un¬ 
certainty in Af(H 2 ) of ~5 x 10^® cm“® ISection 16.21) . 

6.2. Uncertainty in Af(H 2 ) 

As iLee et all (I2012D did, we perform a series of Monte 
Carlo simulations to estimate the uncertainty in Af(H 2 ). 
In these simulations, we assess the errors in Af(HI) and 
Ay, and propagate them together. 

For the uncertainty in Af(HI), we combine two terms in 
quadrature: the error from using a fixed velocity width 
(cHLi), and the error from the correction for high optical 
depth (crHi, 2 )- To estimate cthi.i, we produce 1000 A^(HI) 


images using 1000 velocity widths randomly drawn from 
a Gaussian distribution that peaks at 20 km s“^ with Icr 
of 4 km s“^. The standard deviation of the simulated 
A^(HI) is t hen c omputed for cthi,i. This cthi,! is what 
iLee et aP (j2012l) used as their uncertainty in A^(HI). 
Similarly, to derive crHi, 2 , we generate 1000 A^(HI) im- 
ages by applying the correction to the N (HI) image from 
iLee et all ( 201^ using 1000 combinations of a and b in 
Equation ([6]). These a and b values are again drawn 
from Gaussian distributions whose peaks and widths cor¬ 
respond to the fitted a and b values and their uncertain¬ 
ties. We find that the median Icr of A^(HI) in our study 
is ~8.1 X 10^® cm“®. 

For the uncertain ty in Ay , we repeat the exercise done 
bv iLee et al.l (l2012f ): deriving a number of Ay images by 
changing input conditions (Icr noises of the IRIS 60/100 
fim images, j3, zero point calibration for tiqo), and esti¬ 
mating the minimum and maximum Ay values for each 
pixel. In this exercise, we find that the contribution from 
VSGs to Jeo varies from 76% to 88%, while X varies from 
655 to 855. 

Finally, we propagate the uncertainty in N (HI) and the 
minimum/maximum Ay values through a Monte Carlo 
simulation in order to produce 1000 A^(H 2 ) images. The 
distribution of the simulated A^(H 2 ) is then used to esti¬ 
mate the uncertainty in A^(H 2 ) on a pixel-by-pixel basis. 
We find that the median Icr of A^(H 2 ) in our study is 
^4.8 X 10^® cm“®. 
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Fig. 11.— (left) Rh 2 versus Ejji + Eh 2 for IC348 and BIE (this study). All finite pixels in the rectangular boxes shown in Figure|9]are 
used for p lotti ng. The median 3cr values of i?H 2 and Shi + Sh 2 are indicated as the black dashed lines, while the best-fit curves d etermined 
in Sec tion 16.41 are overlaid in red. The best-fit parameters are shown in the bottom right corner of each plot, (right) Plots from ILee et ahl 


6.3. i?H 2 versus Shi + Sh 2 and Shi versus Shi + Sh 2 

From the rederived iV(HI) and iV(H 2 ) images, we esti¬ 
mate Shi and Sh 2 by 


Shi (Mq pc 


A^(HI) (cm 
1.25 X 1020 


Sh 2 (Mq pc 


iV(H 2 ) (cm 2) 
6.25 X 1019 


(16) 


We find that Shi varies by only a factor of ^^2.6 from 
^4.8 Mq pc“^ to ^12.7 Mq pc“^. On the other hand, 
Sh 2 ranges from —2.4 Mq pc“^ to 81.8 Mq pc“^, al¬ 
though ^98% of the pixels have Sh 2 < 15 Mq pc“^. As 
a result, Shi + Sh 2 bas a small dynamic range of ~8-30 
Mq pc“^ across most of the cloud. 

To compare with the KMT09 predictions aiming at re¬ 
visiting the HI saturation in Perseus, we focus on the in¬ 
dividual dark (B5, BIE, and Bl) and star-forming (IC348 
and NGC1333) regions. The boundaries of each region 
w ere determ ined based on the emission (Section 5 

of iLee et al.l 2012 for details), and are shown as the red 
rectangular boxes in Figure [9l Using all finite pixels in 
the rectangular boxes, we plot Shi and i?H 2 as a function 
of Shi + Sh 2 for each region, and present the results for 
IC348 and BIE in Figures fltH and F T1 Additionally, we 
show the same plots from iLee et al. (l2012l) fo r comp ari- 


Lee et al. 


son. Note that both this study and 
elude negative Sh 2 values by using all finite pixels in the 


3r comp a 

](l20il) 


rectangular boxes. Almost all ('^90%) of these negative 
Sh 2 values fluctuate around zero within uncertainties. 


6.4. Comparison to the KMT09 Predictions 

As in iLee et al.l (|2012ll . the following KMT09 predic¬ 
tions are used to fit the observed i?H 2 vs Shi + Sh 2 pro¬ 
files: 


4tc 

Rn2 = Trr 


1-b 


0.8V'(/)mol 


3'0 L 4tc -I- 3((()i„oi - 1)V'. 


-1 


where 


3 /S 


(ip<^d 


and 


i) = x 

X = 2.3 


4 \ /fH 
2.5+ x 


2.5+ xe’ 
1-F3.1Z'°'365 

(/'CNM 


(17) 

(18) 

(19) 

( 20 ) 


Here is the dust optical depth a spherical cloud would 
have if its HI and H 2 are uniformly mixed, and x is the 
ratio of the rate at which LW photons are absorbed by 
dust grains to the rate at which they are absorbed by H 2 . 
In addition, Scomp is the total gas column density, ad is 
the dust absorption cross section per hydrogen nucleus in 
the LW band, p-n is the mean mass per hydrogen nucleus, 
Z' is the metallicity normalized to the value in the solar 
neighborhood, ^cnm is the ratio of the CNM density to 
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TABLE 1 

Fitting Results for -Rh2 versus Shi + Sh2 


Region 

Best-fit 0CNM 

B5 

8.75 ± 1.35 

IC348 

7.40 ± 0.94 

BIE 

7.28 ± 1.13 

B1 

6.93 ± 1.00 

NGC1333 

5.28 ±0.85 


Note. The uncertainty in 0CNM is estimated from the distribution 
of the simulated 1000 flicNM values. 

the minimum CNM density at which the CNM can be in 
pressure balance with the WNM, and finally (j)mo\ is the 
ratio of the H 2 density to the CNM density. We refer to 
Section 6 of iLee et al.l (|2012ll for a detailed summary of 
the KMT09 model. 

As Equations (17)-(20) suggest, i?H 2 in the KMT09 
model is simply a function of total gas column density, 
metallicity, (/)cnm, and (pmoi, and is independent of the 
strength of the radi ation field in wh ich the cloud is em¬ 
bedded. Following iLee et ^ (|2012h . we adopt Z' = 1 
and (j)nio\ = 10 (fiducial value by KMT09)^^, and con¬ 
strain ^icNM by finding the best-fit model for the ob¬ 
served i?H 2 vs Ehi + Eh 2 • For this purpose, we perform 
Monte Carlo simulations where the uncertainties in i?H 2 
and Ehi + Sh 2 are taken into account for model fitting. 
In these simulations, we add random offsets to i?H 2 and 
Ehi + Eh 2 based on their uncertainties, and determine 
the best-fit curve by setting Z' = 1 and (j)mo\ = 10 and 
finding (/>cnm that results in the minimum sum of squared 
residuals. We repeat this process 1000 times, and esti¬ 
mate the best-fit (/)cnm by calculating the median ^^cnm 
among the simulated (^cnm- The derived (()cnm for each 
region is summarized in Table [l] and the best-fit curves 
are shown in red in Figures [TUI and [TTl 

In the KMT09 model, x measures the relative impor¬ 
tance of dust shielding and H 2 self-shielding, and is pre¬ 
dicted to be ~1 for a wide range of galactic environments. 
In this case, a certain amount of Ehi is required to shield 
H 2 againt photodissociation, and H 2 forms out of HI once 
this minimum shielding column density is obtained. The 
KMT09 model predicts the minimum shielding column 
density of ^10 Mq pc“^ for solar metallicity, and this 
is indeed consistent with what we observe in Perseus: 
Ehi saturates at ^7-9 Mq pc“^ for all five regions. The 
level of the HI saturation changes between the regions 
though, from ~7 Mq pc“^ for B5 to ^ 9 Mm pc“^ for 
NGC1333. In comparison with iLee et al.l (l20l2jl . our Ehi 
values are slightly higher due to the correction for high 
optical depth. This correction brings a closer agreement 
with the KMT09 model. 

The excellent agreement with the KMT09 model is also 
evident from i?H 2 vs Ehi + Eh 2 in Figure [TTJ For all five 
regions, we find that i?H 2 steeply rises at small Ehi+Eh 2 , 
turns over at i?H 2 ~ 1, and then slowly increases at large 
Ehi + Eh 2 - In fact, this common trend on a log-linear 
scale is entirely driven by the almost constant Ehi, and 
therefore it is not surprising to find such a good agree¬ 
ment with the KMT09 model where the Ehi saturation is 
predicted. By fitting the KMT09 predictions to the ob- 

We note that Rh 2 is not sensitive to (jimol- For example, with 
our median (pCNM value of ~7 (Table 1), i?H 2 S't ^HI + Sh 2 = 100 
Mq pc“^ varies by only a factor of ~1.1 for </>n,ol = 10-50. 


served i?H 2 vs Ehi + Eh 2 profiles, we constrain ifiqN M ~ 
5-9, which is consistent with what ILee et al.l (l20I2f) esti¬ 
mated within uncertainties. In the KMT09 model, ^cnm 
determines the CNM density by ncNM = '('CNMUmin 
where is the minimum CNM density at which the 
CNM can be in pressure balance with the WNM, and 
</i’CNM = 5-9 translates into Tcnm = 60-80 K (Equation 
19 of KMT09). This Tcnm range is consistent with what 
we observationally constrained for the HI environment of 
Perseus via the HI absorption measurements (Paper I): 
Tcnm mostly ranges for ^10-200 K, and its distribution 
peaks at ^50 K. Finally, we note that ^!)cnm systemati¬ 
cally decreases toward the southwest of Perseus, reflect¬ 
ing the observed region-to-region variations in Ehi- The 
difference in 0 cnm, however, is not significant, and this 
suggests similar y values for all dark and star-forming 
regions (Equation 20). We indeed find y = ~I.I-I.8, 
in agreement with the KMT09 prediction of comparable 
dust-shielding and H 2 self-shielding for H 2 formation. 

7. OPTICALLY THICK HI: ALTERNATIVE TO THE 
“CO-DARK” GAS? 

Recently, Fukui et al. (2014,2015) suggested that the 
“CO-dark” gas in the Galaxy, referring to the interstellar 
gas undetecta ble either in th e 2I-cm HI and 2.6-mm CO 
emission (e.g.. iBolatto et aH 2013), could be dominated 
by the optically thick HI. As our HI absorption mea¬ 
surements provide an independent estimate of the opti¬ 
cal depth, we can test the validity of their claim against 
Perseus on sub-pc scales. To do so, we utilize our old 
and new Af(H 2 ) images (before and after the correction 
for high optical depth), as well as the CfA Ico data. 

First, we examine how the “CO-dark” gas and the op¬ 
tically thick HI gas are spatially distributed. In order to 
identify the “CO-dark” gas, we use the 3tT contours of the 
ol d A^(H 2 ) and CfA Ico images, following the definition 
by ILee et all (I20I2I1 . These contours are overlaid on our 
correction factor image (Figure^]), and show that the rel¬ 
ative distribution of H 2 (or simply the gas not probed by 
the HI emission) and CO changes across the cloud. For 
example, H 2 and CO trace each other in the southwest, 
while H 2 is more extended than CO elsewhere. ILee et aH 
(I 2 OI 2 II compared H 2 and CO radial profiles for Perseus, 
and estimated that H 2 is on average ^1.4 times more ex¬ 
tended than CO, suggesting a substantial amount of the 
“CO-dark” gas. We then find that the distributions of 
the “CO-dark” gas (traced by the difference between the 
H 2 and CO 3a contours) and the optically thick HI gas 
(traced by the correction factor) generally disagree with 
each other. For example, the region around B5 where 
the discrepancy between the H 2 and CO distributions is 
greatest shows relatively low correction factors. On the 
other hand, the regions with high correction factors at 
(R.A.,decl.) ~ (3‘’33”,-b34°30') and - (3’^23“,-b29°) do 
not have a large amount of the “CO-dark” gas. 

While our previous evaluation was based on the visual 
examination of the 3tT contours, here we more rigorously 
investigate whether or not the optically thick HI gas can 
explain the “CO-dark” gas by comparing the old and new 
A^(H 2 ) images. We first smooth the old and new A^(H 2 ) 
images to 8.4' resolution, as well as their uncertainties, to 
match the resolution of the CfA Ico image. We then find 
the “CO-dark” gas from the smoothed A^(H 2 ) images by 
utilizing the 3cr contours of A^(H 2 ) and Ico- Essentially, 
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logio 2N(Hg) (cm'^) 

Fig. 12.— Histograms of the “CO-dark” gas column density. 
The “CO-dark” gas from the old A^(H 2 ) image is in black, while 
that from the new A^(H 2 ) image is in gray. Both histograms are 
constructed using the data smoothed to 8.4' resolution. 

a pixel is classified as the “CO-dark” gas when 7V(H2) 
is above the 3cr level, but Ico is less than the 3 ct noise. 
For the selected pixels, we calculate the “CO-dark” gas 
column density, 2 iV(H 2 ), and present two histograms in 
Figure IT^ The “CO-dark” gas from the old A^(H 2 ) is in 
black, while that from the new A^(H 2 ) is in gray. Fig¬ 
ure M shows that the two histograms are comparable 
regarding their minimum, maximum, and median val¬ 
ues (different by less than a factor of 2), although the 
gray histogram has a smaller number of pixels due to 
the larger uncertainty of the new iV(H 2 ) image. Given 
that the optically thick HI gas was already taken into 
consideration in the derivation of the new iV(H 2 ) image, 
the comparable black and gray histograms suggest that 
the increased column density due to the optically thick 
HI gas is small (up to ~2 x 10^° cm“^), and the “CO- 
dark” gas still exists in Perseus. In terms of mass, the 
additional contribution from the optically thick HI only 
accounts for ~20% of the observed “CO-dark” gas. 

While our results are in contrast with Fukui et al. 
(2014,2015) who found that the optically thick HI adds 
the column density of ~10^°-10^^ cm“^ and possibly ex¬ 
plains the “CO-dark” gas in the Galaxy, there are multi¬ 
ple factors that could affect the comparison, e.g., spatial 
coverage (^500 deg^ for Perseus vs whole Galactic sky at 
|6| > 15°) and method for deriving A^(H 2 ) (IRIS/2MASS 
vs Planck). In particular, we note that this study and 
Fukui et al. (2014,2015) probe very different scales: our 
results are based on pencil-beam HI absorption measure¬ 
ments on 3.5' scales, while Fukui et al. (2014,2015) esti¬ 
mated the correction factor on 33' scales. If the GNM is 
highly structured with a low filling factor, this could af¬ 
fect the estimate of the correction factor in both studies. 
In the future, it will be important to compare the results 
from Fukui et al. (2014,2015) with a large sample of 
HI absorptio n measurements as w ell a s numerical simu¬ 
lations ('e.g.. lAudit fc Hennebelld 2005: lKim et al.l 2014) 
to investigate how the derivation of the correction factor 
depends on different methodologies and GNM properties. 

8. SUMMARY 

In this paper, we investigate the impact of high optical 
depth on the HI column density distribution across the 
Perseus molecular cloud. We use Arecibo HI emission 
and absorption measurements obtained toward 26 back¬ 


ground sources (Paper I) in order to derive the properties 
of GNM and WNM components along each line of sight 
via the Gaussian decomposition approach (HT03a). The 
derived properties are then used to estimate the correc¬ 
tion factor for high optical depth, and the correction is 
applied to the HI column density image computed in the 
optically thin approximation. 

To revi si t the HI saturation in Perseus observed by 
ILee et al.l (IMll . we rederive the H 2 colum n density 
image by adopting the same methodology as iLee et aH 
(I2012I1 . but using the corrected HI column density im¬ 
age. The final HI and H 2 column density images at ^0.4 
pc resolution are then compared with the KMT09 predic¬ 
tions. Finally, we investigate if the observed “GO-dark” 
gas in Perseus is dominated by the optically thick HI gas. 

We summarize our main results as follows. 

1. We estimate the correction factor for high optical 
depth (/), which is defined as the ratio of the true 
total HI column density (A^tot) to the HI column 
density derived in the optically thin approxima¬ 
tion (A"exp), and express it as a function of Ne^p- / 
= logio(A^exp/102°)(0.32 ± 0.06) -b (0.81 ± 0.05). 
We use two different methods, Gaussian decompo¬ 
sition and isothermal approximation methods, and 
find that they are consistent within uncertainties. 
This is likely due to the relatively low optical depth 
and insignificant contribution from the diffuse ra¬ 
dio continuum emission for Perseus. 

2. We estimate that the correction factor in/around 
Perseus is small (up to ~1.2), and the total HI mass 
increases by only ~10% from ~2.3 x 10'^ Mq to 
~2.5 X 10^ Mq due to the inclusion of the optically 
thick HI gas. 

3. The H 2 column density image rederived using the 
corrected HI column density image is comparable 
to the original one by ILee et aP (|2012[ 1 , confirming 
the minor correction for high optical depth. 

4. For individual dark and star-forming regions in 
Perseus (B5, BIE, Bl, IC348, and NGC1333), the 
HI surface density is relatively uniform with ~7- 

9 Mpi pc ~^. _ This is slightly higher than what 

ILee et all (|2ni2h found due to the correction for 
high optical depth. The correction brings a closer 
agreement with the KMT09 model where the min¬ 
imum HI surface density of ^10 Mq pc“^ is pre¬ 
dicted for shielding H 2 against photodissociation in 
the solar metallicity environment. 

5. Driven by the uniform Ehi ~ 7-9 Mq pc“^, i?H 2 vs 
Shi + Sh 2 on a log-linear scale shows remarkably 
consistent results for all dark and star-forming re¬ 
gions: i?H 2 sharply rises at small Ehi + Sh 2 , and 
then gradually increases toward large Ehi + Eh 2 
with the transition at i?H 2 ^ 1- 

6. The mass increase due to the optically thick HI only 
accounts for ^^20% of the observed “CO-dark” gas 
in Perseus. The spatial distributions of the “CO- 
dark” gas and the optically thick HI gas do not 
generally coincide with each other, and the “CO- 
dark” gas still exists even after the optically thick 
HI is considered in the derivation of H 2 . 

Our study is one of the first attempts to ex¬ 
amine the properties of the cold and warm HI 
in molecular cloud environments and their relation 
with the HI-to-H 2 transition. While HI envelopes 
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surrounding molecula r clouds have been f requently 
found (e.^. iKnaod 1 974: IWannier el alJ 1983; 


Elmegreen fc Elmeereenl 1987; Andersson et al.l 1991; 


Rogers et al.l 199 5: IWilliams &: Maddalena 1996; 
Imara fc Blitd 2011; iLee et aL 2012: iLee et al. 2014; 

2014), HI has traditionally been considered 


Motte et al 


less important for the formation and evolution of molec¬ 
ular clouds. The excellent agreement between our ob¬ 
servations of Perseus on sub-pc scales and the KMT09 
model, on the other hand, suggests the significance of HI 
as one of the key ingredients for the Hl-to-H 2 transition 
and consequently for star formation. 

While our data agree with the KMT09 model, further 
theoretical developments are required. Eor example, the 
KMT09 model considers only the CNM as a source of 
shielding against H 2 photodissociaiton. Our HI emission 
and absorption measurements, however, show that the 
CNM and the WNM have comparable column density 
contributions for Perseus (Paper 1). Clearly, the WNM 
needs to be taken into consideration in the models of 
H 2 formation (e.g., Bialy et al. in prep). In addition, 
the KMT09 model ignores the effect of internal radi¬ 
ation field. This is valid for Perseus, as there are no 
early-type stars p roducing a significant amount of disso¬ 
ciating radiation (ILee et al.l 2012). However, the role of 
internal radiation field would be important for massive 
molecular clouds containing a large number of OB stars. 
The discrepancy with the KMT09 model recently found 
for th e W43 molec ular cloud complex in the Galactic 
plane (jMotte et al.l 2014) and the low-metallicity SMC 


(|Weltv et al.l 2012) suggests that additional physical in¬ 
gredients (e.g., better description of H 2 formation and 
photodissociation in the low-metallicity ISM and strong 
shear/turbulence at galactic bars) would be necessary for 
extreme environments. Einally, vertical thermal and dy¬ 
namical equilibrium in a galactic disk is another impor¬ 
tant aspect, as recently explored bv iKim et al.l (120141) . 
Euture comparisons between theoretical models and HI 
emission/absorption observations of molecular clouds in 
a wide range of ISM environments will be important for 
a deep understanding of HI properties and their role in 
the formation and evolution of molecular clouds. 
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APPENDIX 


A. COMPARISON BETWEEN METHOD 1 AND METHOD 2 

In order to investigate why the Gaussian decomposition and isothermal methods result in comparable correction 
factors despite of their apparent differences ('Sections 14.21 and 14.31) . we start with Equation (fT4|l . and compare it with 
Equation ([5]): 

Ntot (cm“^) = Co y /chan(l^)Texp(w)c?U (Al) 


where Cq = 1.823 x 10^® cm ^/(K km s ^). Gombined with Equation (TT^ . Equation (jAip then becomes 

Ntot (cm“^) = Co [ — , ^ - ^ ^ , Tg:x.n{v)dv 


(A2) 


The approximation in Equation (IA2() can be made due to the small Ts{v)/[Ts{v) — Tsky] of ^1.0-1.2 lSection l4Al) . 

While an optical depth spectrum is primarily determined by the GNM, an “expected” emission spectrum has con¬ 
tributions from both the CNM and the WNM. Hence, Texp{v) can be expressed as (Equations 2, 3, and 4 of Paper 

I): 

'Texp('y) = T'b,CNm(t) + 'rB,WNM(l') (A3) 


with 


N-l 

TB,CNMiv) = ^ r,,„(l 
0 

K-1 

rB.WNM(u) = ^ [Ffc + (1 - 

0 


where rB,CNM(i^) is the emission from N CNM components, T'b,wnm('I’) is the emission from K WNM components, 
the components with subscript n refer to the CNM, the components with subscript k refer to the WNM, is the 
optical depth of m CNM component that lies in front of n cloud, Fk is the fraction of k WNM component that lies 
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in front of all CNM clouds, Tq is th e pe ak brightness temperature, vq is the central velocity, and finally 6v is the 1/e 
width of the component. Equation (IA2I) then becomes 


, , (N-l K-l 'j 

Ntot (cm- 2 ) r^Co ^ + (1 - I dv. 

® I 0 0 J 

(A4) 

We examine Equation (IA4|) under two extreme circumstances: (I) F^. = 0 and r <C 1 and (II) E/ = 1 and t <C 1. In 
both cases, the first part of Equation (IA4I) becomes 

A = Co 


Cn 


Co 


t{v) 

l"y 

1 - e-^G) 

It 

t{v) 

1^ 

1 - e-^G) 

It 

t(u) 



dv 


. N-l 


= C'o / ^ Ts^nTn{v)dv 

d 0 

= Co [ Ts,nT0,neAC-o,A/Sv.]Av. 

^ n 


On the other hand, the second part of Equation (IA4I1 becomes 
(I) Fk = 0 (all WNM components are behind the CNM) and t <C 1: 


B = C, 


= Q 


C, 


I I AA) I E I-f’* + (' - I dv 

K-l ^ 


Tjv) 

1-(1-t(u)) 


» K-l 

n 


(II) Fk = 1 (all WNM components are in front of the CNM) and r <C I: 

B = Co 


= Cr 


Cr 


/ I - I E | dv 

K-l ^ 


t{v) 


I - e-^G) 

It 

t{v) 


1 - e-'^C) 

It 

t{v) 



1-(1-t(u)) 


„K-1 

n 


As a result, we have 

Ntot (cm“^) = A + B 


C, 


. (N-l K-l -) 

/ E 7^s.nro.„e-[(—+ E dv, 

I 0 0 J 


(A5) 


(A6) 


(A7) 


(A8) 


which is essentially Equation 

These examples of the two extreme cases suggest that the Gaussian decomposition and isothermal methods would 
result in comparable correction factors (/ = A^tot/A^exp) regardless of Fk when r <C 1 and Tgky Tg. 
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Fig. 13.— / = A^tot/-^exp Q-s a function of logio(A^exp/10^^) for 68 sources in HT03a,b. 11 lines of sight probing the Galactic plane (|6| < 
10°) are excluded from the analysis to make a comparison with our study of Perseus. Both the true and optically thin HI column densities 
{Ntot and A^exp) are calculated using the published full line of sight information, and the linear fit to all 68 data points is overlaid in red: 
/ = logio(7Vexp/lo2°)(0.26 ± 0.08) + (0.91 ± 0.04). 

B. CORRECTION FACTOR IN HEILES & TROLAND (2003A,B) 

To examine the validity of the linear relation we assume for the correction factor and the optically thin HI column 
density (Sections 14.21 and l4..'l|l . we produce the same plot as Figured (left) for HT03a,b. Among the total 79 sources 
in HT03a,b, 11 sources with |6| < 10° are excluded for a fair comparison with our study. The published “expected” 
emission spectra and total HI column densities are used to estimate / and logio(lVexp/10^°) in the same way as we do 
for our 26 sources, and the results are presented in Figure [131 
We find that Figure [T51 is very similar with Figures [3] (left) and |S] (left), suggesting that our assumption of the 
linear relation between / and logio(lVexp/10^°) for Perseus is reasonable. Interestingly, however, a few sources with 
logio(iVexp/10^°) > 1 show some deivation from the linear relation, which requires a further examination. See Section 
14.31 for a more discussion. 
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